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ABSTRACT 

We report the discovery of 1RXS J173006. 4+033813, a polar cataclysmic variable with a period of 
120.21 min. The white dwarf primary has a magnetic field of B = 42+g MG, and the secondary is a 
M3 dwarf. The system shows highly symmetric double peaked photometric modulation in the active 
state as well as in quiescence. These arise from a combination of cyclotron beaming and ellipsoidal 
modulation. The projected orbital velocity of the secondary is K% — 390 ± 4kms _1 . We place an 
upper limit of 830 ± 65 pc on the distance. 

Subject headings: binaries: close — binaries: spectroscopic — novae, cataclysmic variables — stars: 
individiual (1RXS J173006.4+033813) — stars: variables: other 



1. INTRODUCTION 

Cataclysmic Variables (CVs) are close interacting bi- 
nary systems in which a white dwarf (WD) accretes ma- 
terial from a Roche lobe fil ling late-type secondary star 
(|Warnerlll995l: lHellierll2001h . In most non-magnetic CVs 
(B < 10 4 G), the material lost from the secondary does 
not directly fall onto the WD because of its large spe- 
cific orbital momentum: instead, it settles down in an 
accretion disc around the WD. 

The accretion disc is the brightest component of the 
CV due to the large gravitational energy release in vis- 
cous accretion. The disc dominates the emission from 
the WD and donor over a wide wavelength range. 

On the other hand, the accretion geometry in mag- 
netic CVs is strongly influenced by the WD magnetic 
field. Magnetic CVs are broadly divided into two sub- 
classes: Polars and Intermediate Polars (IPs). Polars 
usually show a synchronous or near synchronous rotation 
of WD with the orbital motion of the binary system and 
have high mag netic fields (B > 10 MG) (for a review, see 
ICropper] [l990). In IPs the WD rotation is far from syn- 
chronous and ty pically have magnetic field, B < 10 MG 
(for a review, see lHellierll2002l ). The strong magnetic field 
in polars deflects the accretion material from a ballistic 
trajectory before an accretion disc can form, channeling 
it to the WD magnetic pole(s). The infalling material 
forms a shock near the WD surface, which produces ra- 
diation from X-rays to infrared wavelengths. Electrons 
in the ionized shocked region spiral around the magnetic 
field lines and emit strongly polarized cyclotron radia- 
tion at optical and infrared wavelengths. Polars exhibit 
X-ray on (high) and off (l ow) states more frequ ently than 
the other variety of CVs (|Ramsav et al.ll2004D . 

1RXS J173006.4+033813 (hereafter 1RXS J1730+03) 
is a Galactic source that is highly variable in the op- 
tical and X-ray, exhibiting dramatic outbursts of more 
than 3 magnitudes in optical. It was disco vered by the 
ROSAT satellite during its all-sky survey (jVoges et al.l 

1 Department of Astronomy, California Institute of Technol- 
ogy, Pasadena, CA 91125, USA 

On sabbatical leave from Department of Astronomy and As- 
trophysics, University of Toronto, 50St. George Street, Toronto, 
ON M5S 3H4, Canada. 

3 Space Radiation Laboratory, California Institute of Technol- 
ogy, Pasadena, CA 91125, USA 



1999). iDenisenko et al.1 (|2009T ) . in the course of their in- 
vestigation of poorly studied ROSAT sources, reported 
that USNO-B1.0 object 0936-00303814 which is within 
the 10" (radius) localization of the X-ray source showed 
great variability (AR of up to 3 mag) in archival data 
(Palomar Sky Survey; SkyMorph/NEAT). During cer- 
tain epochs the s ource appears to have be en undetectable 
(rriR > 20 mag). IDenisenko et al.l ()2009l ) undertook ob- 
servations with Kazan State University's 30-cm robotic 
telescope and found variability on rapid timescales of 10 
minutes. 

In this paper, we report the results of our photometric, 
spectroscopic and X-ray follow-up of 1RXS J1730+03. 

2. OBSERVATIONS 
2.1. Optical photometry 

We observed 1RXS J1730 +03 with the Pa lomar 
Robotic 60-inch telescope (P60: ICenko et aLll2006l ) from 
UT 2009 April 17 to UT 2009 June 5, and w ith the Large 
Format Camera (LFC; iSimcoe et al.1 120001 ) at the 5 m 
Hale telescope at Palomar on UT 2009 August 26. Here 
we give details of the photometry. 

We define a photometric epoch as observations from 
a single night when the source could be observed. We 
obtained 28 epochs with the P60, subject to schedul- 
ing and weather constraints. A typical epoch consists of 
consecutive 90-120 s exposures spanning between 30-300 
minutes (Table [TJ. We obtained <?', r', i' photometry on 
the first and third epochs. After the third epoch, we 
continued monitoring the source only in i' band. 

We reduced the raw images using the default P60 im- 
age analysis pipeline. LFC images were reduced in IRAfQ. 
We performed photom etry using the IDI@ DADPHOT pack- 
age ((Landsman 1993]) ■ Fluxes of the target and reference 
stars (Figure [l} were extracted using the APER routine. 
For aperture photometry, the extractio n region was se t 
to one seeing radius, as recommended bylMighedlj (1999). 
The sky background was extracted from an annular re- 
gion 5-15 seeing radii wide. We used flux zero points and 
seeing values output by the P60 analysis pipeline. Mag- 
nitudes for the reference stars were calculated from a few 
images. The magnitude of 1RXS J1730+03 was calcu- 

4 http://iraf.noao.edu/ 

5 http:/ /www. ittvis.com/ProductServices/IDL.aspx 



2 



Bhalerao et. al 



Table 1 

Photometry of 1RXS J1730+03 



Date (UT) 


HJD 


Filter name 3 


Exposure time (sec) 


Magnitude 


Error b 


20090417 


54938.832278 


9' 


120 


20.8 


0.13 


20090417 


54938.835488 


i' 


90 


19.39 


0.05 


20090417 


54938.838696 


g' 


90 


20.75 


0.13 


20090417 


54938.839953 


%' 


90 


19.17 


0.05 


20090417 


54938.841212 


g' 


90 


20.67 


0.13 


20090417 


54938.842469 


i' 


90 


19.1 


0.05 


20090417 


54938.843727 


g' 


90 


20.63 


0.12 


20090826 


55069.691401 


Ifci' 


60 


20.4 


0.05 


20090826 


55069.693218 


Ifci' 


60 


20.04 


0.06 



Note. — This table is available in a machine readable form online. A part of the table is reproduced here for 
demonstrating the form and content of the table. 

a Filters g' , r' and i' denote data acquired at P60 in the respective filters, Ifci' denotes data acquired in the i' band 
with the Large Format Camera at the Palomar 200" Hale telescope 

b Relative photometry error. Values do not include an absolute photometry uncertainty of 0.16 mag in the g' band, 
0.14mag in the r' band and 0.06mag in the i' band. Absolute photometry is derived from default P60 zero point 
calibrations. 



Table 2 

Photometry of reference stars for 1RXS J1730+03 



Identifier 51 


Right Ascension 


Declination 


g 1 magnitude 


r' magnitude 


i' magnitude 


A 


262:30:21.99 


03:38:37.5 


15.861 ± 0.003 


15.600 ± 0.003 


15.341 ± 0.003 


B 


262:30:14.57 


03:37:11.0 


17.372 ± 0.009 


17.153 ± 0.006 


16.872 ± 0.008 


C, 10 


262:32:45.57 


03:37:17.0 


16.879 ± 0.006 


16.691 ± 0.005 


16.423 ± 0.006 


D, 9 


262:32:13.84 


03:38:03.0 


18.161 ± 0.016 


17.929 ± 0.010 


17.660 ± 0.013 


E 


262:31:00.31 


03:38:31.3 


18.518 ± 0.020 


17.462 pm 0.007 


16.811 pm 0.007 


F 


262:30:55.06 


03:37:26.4 


18.530 ± 0.022 


18.158 pm 0.012 


17.816 pm 0.016 


G, 8 


262:31:57.98 


03:38:03.9 


18.621 ± 0.022 


18.259 pm 0.012 


17.905 pm 0.016 


H 


262:30:41.22 


03:38:32.8 


16.770 ± 0.006 


16.081 pm 0.004 


15.688 pm 0.004 


1 


262:29:45.31 


03:38:43.7 


17.661 ± 0.011 


17.302 pm 0.007 


16.969 pm 0.008 


1 


262:32:01.97 


03:40:28.3 






17.059 ± 0.035 


2 


262:31:44.94 


03:39:48.6 






19.855 ± 0.043 


3 


262:32:31.60 


03:39:32.4 






19.148 ± 0.039 


1 


262:32:33.94 


03:38:59.4 






17.622 ± 0.033 


5 


262:32:11.15 


03:38:39.0 






20.777 ± 0.049 


6 


262:31:35.94 


03:38:29.5 






19.843 ± 0.026 


7 


262:31:07.21 


03:37:57.9 






19.057 ± 0.026 


11 


262:31:55.09 


03:36:45.4 






17.471 ± 0.030 


12 


262:31:40.33 


03:36:34.5 






16.339 ± 0.044 


13 


262:33:19.44 


03:36:21.2 






18.006 ± 0.030 


14 


262:31:45.44 


03:36:12.4 






18.085 ± 0.048 


15 


262:32:18.64 


03:35:48.1 






17.392 ± 0.055 



Note. — This table is available in a machine readable form online. A part of the table is reproduced here for demonstrating the 
form and content of the table. 

a The letters A-I deno te stars used in photometry of P60 data, numbers 1-15 denote reference stars used in photometry of LFC images 
(Figure [Q Section [2~H . 

b Relative photometry error. Values do not include an absolute photometry uncertainty of 0.16 mag in the g' band, 0.14mag in the r' 
band and 0.064 mag in the i' band. Absolute photometry is derived using default P60 zero point calibrations. 



latcd relative to the mean magnitude of a 9 reference 
stars for LFC images, and 15 reference stars for P60 im- 
ages (Table [2]). The LFC images (Figure [J) resolve out a 
faint nearby star (m^ = 20.8), 3". 4 from the target. The 
median seeing in P60 data is 2".l (Gaussian FWHM): so 
there is a slight contribution from the flux of this star 
to photometry of 1RXS J1730+03. We do not correct 
for this contamination. The statistical uncertainty in 
magnitudes is ~ 0.2 mag for P60 and ~ 0.05 mag for 
LFC, and the systematic uncertainty is 0.16 mag in the 
g' band, 0.14mag in the r' band and 0.06 mag in the i' 
band. 

The resultant lightcurves are shown in Figures [21 [3] & 
01 Table Q] provides the photometry. 



2.2. Spectroscopy 

We obtained optical and near-infrared spectra of 
1RXS J1730+03 at various stages after outburst (Fig- 
ure [5]). The first optical spectra were taken 13 days after 
the first photometric epoch. We used the Low Reso- 
lution Imagin g Spectrograph on the 10 m Keck-I tele- 
scop e (LRIS;lOke et al.lll995|), with upgraded blue cam- 
era (jMcCarthv et alj|1998t iSteidel et al.ll2004D . covering 
a wavelength range from 3,200 A -9,200 A. We acquired 
more optical data 34 days after outburst, with the Double 
Beam Spectr ograph on the 5 m Hale telescope at Palo- 
mar (DBSP; lOke k Gunnl [19821) . We took 5 exposures 
spanning one complete photometric period, covering the 
3,500 A -10,000 A wavelength range. We took late time 
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Figure 1. Finder chart for 1RXS J1730+03 (a = 17 h 30 m 06M9, 8 = +03°38'18".8). This i' band image was acquired with the Large 
Format Camera (LFC) at the 5 m Hale telescope at Palomar. Stars numbered 1-15 in red are used for relative photometry in LFC data. 
Stars labeled A— I in blue are used for relative photometry in P60 data (Scction l2.1[ l, The green circle shows the 5" extraction region used 
for calculating UVOT fluxes. The circle includes a contaminator, 3". 4 to the South- West of the target. 



spectra covering just over one photometric period for the 
quiescent source with the upgraded LRIS0. At this epoch, 
we aligned the slit at a position angle of 45 degrees to 
cover both the target and the contaminator, 3". 4 to its 
South West (Figured]). We also obtained low resolution 
J-band spectra with the Near InfraRed Spectrograph 
on th e 10 m Keck-II telescope (NIRSPEC: lMcLean et all 
1998). 12 spectra of 5 minutes each were acquired, cover- 
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ing the wavelength region from 11,500 A- 13,700 A. For 
details of the observing set up, see the notes to Table [5j 
We analyzed the spectra using IRAF and MIDASQ and 
flux calibrated them using appropriate standards. Wave- 
length solutions were obtained using arc lamps and with 
offsets determined from sky emission lines. Figure [S] 
shows an optical spectrum from each epoch, while the 
IR spectrum is shown in Figure [71 

Munich Image Data Analysis System; 

/ww. eso.org/sci/data-proccssing/softwarc/esomidas/ 
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Table 3 

Observation Log for Swift ToO Observations of 1RXS J1730+03. 



Obs ID 


Start Date & Time 


Stop Time 


Exposure 


Filter 


Wavelength a 


Magnitude 


Flux 


Flux 








(s) 




(A) 






(A*Jy) 


00035571001 


2006 Feb 9 16:56:43 


18:42:00 


1106 


UVM2 


2231 


17.9±0.1 


31±2 b 


51±3 










XRT 






0.02 c 


0.06 (5keV) 


00031408001 


2009 May 3 18:28:56 


19:02:12 


1964 


UVM2 


2231 


20.44±0.21 


3.04±0.59 b 


5.05±0.99 


00031408002 


2009 May 4 18:34:04 


22:12:54 


4896 


UVW1 


2634 


20.14±0.10 


3.47±0.32 b 


8.04±0.75 


00031408003 


2009 May 6 03:02:21 


15:36:33 


5621 


UVW2 


2030 


21.48±0.17 


1.38±0.22 b 


1.93±0.30 








XRT 






<0.002 c 


<0.006 (5keV) 



a Effective wavelength for each filter for a Vega-like spectrum 



b Flux in the units of 10 _1 
c Counts s" 1 in 0.5-10 keV. 
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Figure 2. P60 photometry of 1RXS J1730+03. Top panel: Epoch 
1 photometry in g' and i' bands. Bottom panel: Epoch 3 photom- 
etry in r' and i' bands. 
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The second LRIS epoch had variable sky conditions. 
Here, we extracted spectra of the aforementioned con- 
taminator. This object is also a M-dwarf, hence both 
target and contaminator spectra will be similarly affected 
by the atmosphere. We estimate the i' magnitude of the 
contaminator in each spectrum, and compare it to the 
the value measured from the LFC images to estimate 
and correct for the extinction by clouds. 

2.3. X-ray and UV observations 

We observed the 1RXS J1730+03 with the X-ray tele- 
scope (XRT) and the UV-Optical Telescope (UVOT) 
onboard the Swift X-ray satellite ()Gehrels et al.l 120041) 
during UT 2009 May 3-6 for a total of about 12.5 ks. 
The level-two event data was processed using Swift data 
analysis threads for the XRT (Photon counting mode; 
PC) using the HEASARC FT00L^| software package. The 
source was not detected in the X-ray band. 

We follow the procedures outlined by iPoole et al.1 
(2008) for analyzing the UVOT data. The measured 
fluxes are given in Table [3J The contaminator is within 
the recommended 5" extraction radius. Hence, flux mea- 
s urements are upper lim its . 

IShevchuk et aLl (l2009l) had observed 1RXS J 1730+03 
on UT 2006 February 9 with the Swift satellite as a part 
of investigations of unidentified ROSAT sources. They 
detected the source with a count rate of 0.02 counts s . 
The best-fit power law has a photon index T — 1.8 ± 
0.5 and a 0.5-10keV flux of 1.2 x 10~ 12 ergcnT 2 s" 1 . 
After converting to the ROSAT bandpass assuming the 
XRT model parameters, this value is approximately a 
factor of two lower than the archival ROSAT flux. They 
report a much higher UV flux in their observations, which 
suggests that the source was in an active state during 
their observations. 

The column density inferred from the XRT data is low, 
JV H ~ 7 x 10 20 cm" 2 . From ROSAT data, this column 
densi ty correspo nds to Ay = 0.39, which gives Atjym2 = 
0.84 (jCoxl I200QO . For comparison. iSchlegel et all (1998) 



give the Galactic dust extinction towards this direction 
(I = 26. °7, b = 19°. 7) to be E(B - V) = 0.141 mag 
(Ay = 0.44), corresponding to a column density of about 
8xl0 20 cnr 2 . 

3. NATURE OF THE COMPONENTS 

The optical spectra (Figure [6]) show rising flux towards 
the red and blue ends of the spectrum: indicative of a 

8 http://heasarc.gsfc.nasa.gov/ftools/; Blackburn (1995) 
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Figure 3. i' band P60 photometry of 1RXS J1730+03, for Epochs 1 to 28. The source was in outburst in Epoch 1 and sl owly faded into 
quiescence. The error bars denote relative photometry errors and do not include an absolute zero point error of 0.064 mag (§ 12.11) , The solid 
blue line denotes the best fit sinusoid to the photometric data. The dashed green line is added to show the relative phase of the source in 
outburst, these epochs were not used for fitting the lightcurve. 



hot (blue) and cool (red) component. The red part of 
the spectrum shows clear molecular features, character- 
istic of late type stars. The blue component is devoid of 
any prominent absorption/emission features. From the 
overall spectral shape we infer that 1RXS J1730+03 is a 
CV. 

3.1. Red component 

The red component of 1RXS J 1730+03 is typical of 
a late type star. In Figure [8j we compare the red 
side spectrum of 1RXS J1730+03 with several M-dwarfs. 



From the shape of the TiO bands at 7053 A -7861 A, 
we infer that the spectral type to be M3±l. This is 
consis tent with the relati vely featureless J band spec- 
trum (Mc Lean et al.ll2003T) . The s pectral typ e indicates 
an effective temperature of 3400 K (Cox 2000). The pres- 
ence of a sodium doublet at 8183/8195A implies a lumi- 
nosity class V. 

We also f it the spectrum with model atmospheres cal- 
culated by iMunari et al.l |2005). For late type stars, 
these models are calculated in steps of AT = 500 K, 
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Table 4 

Locations of Cyclotron Harmonics. 



Harmonic 


Measured 


Measured 


Interred 


number 


Wavelength 


Frequency 


Wavelength a 




(A) 


(Hz) 


(A) 


7 


3540 


8.4 x 10 14 


3664 


6 


4440 


6.7 x 10 14 


4275 


5 


5180 


5.7 x 10 14 


5130 


4 


6770 


4.4 x 10 14 


6413 


3 


8225 


3.6 x 10 14 


8551 


2 






12826 


1 






25653 



Calculated from the best fit magnetic field strength, B = 42 MG 



Figure 4. i' band LFC photometry of 1RXS J1730+03 in quies- 
cence, 131 days after we first saw high variability. The variations 
have the same phase and period as observed during outburst. See 
Figure l3l for details. 



A log g = 0.5 and A[M/H] = 0.5. We use model at- 
mospheres with no rotational velocity (V ro t = Okms -1 ) 
and convolve them with a kernel modeled on the seeing, 
slit size and pixel size. We ignore the regions contami- 
nated by the telluric A and B bands (7615 A and 6875 A). 
The unknown contribution from the white dwarf was fit 
as a low order polynomial. We correct for extinction 
using Ay = 0.39 from X-ray data (Section l2\3|) . To mea- 
sure log g, we use the spectrum in the 8,000 A -8,700 A 
region, which is expected to have fairly little contami- 
nation from the blue component. This region includes 
the Ca II lines at 8498, 8542 A and the Na I doublet, 
which are sensitive to log g. We then fit the spectra 
m the 6,700 A -8,700 A range to determine the tempera- 
ture and metallicity. The best fit model has T = 3500 K, 
log g = 5.0 and solar metallicity, consistent with our de- 
t ermination of the s pectral type. 

iKolb et ail (|200lD state that unevolved donors in CVs 
follow the spectral type-mass relation of the zero age 
main sequence, as the effects of thermal disequilibrium 
on the secondary spectral type are negligible. For a M3 
star, this yields a mass of 0.38 Mq. As the secondaries 
evolve, the spectral type is no longer a good indicator of 
the mass and gives only an upper limit on the mass. The 
lower limit can simply be assumed to be the Hydrogen- 
burning limit of 0.08 Mq. 

3.2. Blue component 

The blue component of 1RXS J1730+03 is consistent 
with a highly magnetic white dwarf. The blue spec- 
trum is suggestive of a hot object, but does not show 
any prominent absorption/emission features. Ha is seen 
in emission, but other Balmer features are not detected. 
The spectrum (Figure [6| shows cyclotron humps, sug- 
gesting the presence of a strong magnetic field. The po- 
lar nature of the object is supported by the absence of 
an accretion disc, and the transition from an active state 
to an off state in X-rays (Section[TJ Section l2"U)) . 

For analyzing the WD spectrum, we subtracted a 
scaled spectrum of the M-dwarf GL694 from the com- 
posite spectrum of 1RXS J1730+03. The resultant spec- 
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Figure 5. Photometric evolution of 1RXS J1730+03 in i' band. 
Data for each epoch was fit with a sinusoid with the same period 
for all epochs. In contrast to Figures [3] & H] the phase was allowed 
to vary independently for each epoch. Top: average i' magnitude 
for each P60 observation epoch, as a function of time. Blue arrows 
mark spectroscopy epochs. Bottom: semi-amplitude of sinusoidal 
variations as a function of the corresponding mean i 1 magnitudes 
for each epoch. 
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are diagonally hatched. The first spectrum was obtained 13 days after we measured high variability, the second after 34 days, and the 
lowermost spectrum after 60 days. The spectra show broad Balmer features which evolve with time. The M-dwarf features (TiO, Na I) are 
clearly seen at all times. 




Figure 7. Keck NIRSPEC spectrum of 1RXS J1730+03. The spectrum is fiat and nearly featureless, as expected for early M stars. 



trum (Figure [9]) clearly shows cyclotron harmonics. The 
hump seen in the J-band spectrum (Figure [7]) is also 
inferred to be a cyclotron harmonic. A detailed mod- 
eling of the magnetic field is beyond the scope of this 
work, but we use a simple model to estimate the mag- 
netic field. We fit the cyclotron humps with Gaus- 
sians and measure the central wavelengths (Table 2]). 
We then fit these as a series of harmonics, and infer 
that the cyclotron frequency is v cyc = 1.17 x 10 14 Hz. 
The magnetic field in the emission region is given by, 
B = (^cyc/2.8 x 10 14 Hz) • 10 8 G = 42 x 10 6 G. 

As a conservative error estimate, we consider the worst 
case scenario where our identified locations for the cy- 
clotron humps arc off by half the spacing between con- 
secutive cyclotron harmonics. Using this, we estimate 
the errors on the magnetic field: B = 42+g MG. 

4. SYSTEM PARAMETERS 
4.1. Orbit 

We use the best fit iMunari et~a l. (2005) model atmo- 
sphere to measure radial velocities of the M-dwarf. We 
vary the radial velocity of the model, and minimize the 
X 2 over the 6,500 A -8,700 A spectral region, excluding 



the telluric O2 bands. After a first iteration, the spectra 
are re- fit to account for motion of the M-dwarf during the 
integration time. For the 12 spectra taken at the second 
LRIS epoch, we also measure the radial velocity for the 
contaminator star on the slit, and find it to be constant. 
This serves as a useful test for our radial velocity mea- 
surement procedure. The barycentric corrected velocities 
are given in Table [5] 

We fit a circular orbit (v 2 = 72+^2 sm([2n(t — t )]/P)) 
to the measured velocities. We define the superior con- 
junctioifl of the WD as phase 0. 

The 2009 June 16 spectroscopic data (Table [5]) give 
an orbital period P C st = 123 ± 3min. We then use the 
photometric variability (Section l4.2|) to determine an ac- 
curate period in this range. Next, we refine the solution 
with velocity measurements from the other two spec- 
troscopic epochs. The best-fit solution gives a period 
P = 120.2090 ± 0.0013 min and K 2 = 390 ± 4kms~ 1 
(Table Figure [TO]). 

Figure [TT] shows sections of the last epoch spectra 

9 When the WD is furthest from the observer along the line of 
sight. 
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Figure 8. Comparison of red part of 1RXS J1730+03 spectrum 
with M-dwarf spectra. Some spectra are not corrected for Telluric 
absorption, the affected regions are diagonally hatched. Prominent 
bands (TiO, CaH) and lines (Ha, Na I) are marked. Comparing 
the shape of the TiO bands at 7053-7861 A and the shape of the 
continuum red wards of 8200 A, we infer that spectral type of the 
red component to be M3. The presence of a sodium doublet at 
8183/8195 A implies a luminosity class V. 



Table 5 

Radial velocity of the M-dwarf. 



Heliocentric JD Exposure time Barycentri Instrument 

Radial Velocity 





(s) 


(kms 1 ) 




2454952.08721 


1020 


—398 ± 34 


LRIS a 


2454952.10067 


1020 


— 195 ± 22 


LRIS a 


2454952.11179 


600 


109 ± 16 


LRIS a 


2454972.92923 


600 


55 ± 92 


DBSP b 


2454972.93818 


900 


-258 ± 112 


DBSP b 


2454972.94896 


900 


-324 ± 69 


DBSP b 


2454972.95956 


900 


-430 ± 105 


DBSP b 


2454972.96875 


600 


-42 ± 150 


DBSP b 


2454998.99566 


210 


-421 ± 11 


LRIS C 


2454998.99969 


300 


-438 ± 9 


LRIS C 


2454999.00418 


300 


-417 ± 9 


LRIS C 


2454999.00850 


300 


-361 ± 9 


LRIS C 


2454999.01283 


300 


-266 ± 10 


LRIS C 


2454999.01715 


300 


-163 ±12 


LRIS C 


2454999.02147 


300 


-32 ± 12 


LRIS C 


2454999.02580 


300 


98 ± 15 


LRIS C 


2454999.03012 


300 


186 ± 10 


LRIS C 


2454999.03448 


300 


290 ± 17 


LRIS C 


2454999.03880 


300 


336 ± 13 


LRIS C 


2454999.04317 


300 


344 ± 13 


LRIS C 



a Low Resol ution Ima g ing S pectrograph on the 10 m Keck- 
I t elescope (Okc ct al. 1995) , with up graded blue cam- 
era (McCartK^t_al. 1998; Steidcl ct al. 200|). Settings: Blue side: 
3,200 A -5,760 A, grism with 400 grooves/mm, blaze 3,400 A, dis- 
persion l.OQApixel- 1 , R -700. Red side: 5450A-9250A, grat- 
ing with 400 grooves/mm, blaze 8,500 A, dispersion 1.16 A pixel -1 , 
R ~ 1,600. Dichroic: 5,600 A, slit: 1".0. 

b Double Spectrograph on the 5 m Hale telescope at Palo- 
mar l|Oke fc Gunnl ll"98"2h . Settings: Blue side: 3,270 A- 
5700 A, grating with 600 lines/mm, blaze 4,000 A, dispersion 
1.08 A pixel -1 , R -1,300. Red side: 5,300 A-10,200 A, grat- 
ing with 158 grooves/mm, blaze 7,500 A, dispersion 4.9Apixel — 1 , 
R ~600. Dichroic: D55 (5500A), slit: 1".5. 

c Low Resolution Imaging Spectrograph on the 10 m 
Keck-I telescope, with upgraded blue and red cameras: 
|http:/ /www2. keck, hawaii.edu/inst/lris/lris-red-upgrade-notes . html 
Settings: Blue side: 3,300 A — 5, 700 A, grism with 400 grooves/mm, 
blaze 3,400 A, dispersion 1.09 Apixel -1 , R ~ 700. Red side: 
6,760A-8,800A, grating with 831 grooves/mm, blaze 8200 A, 
dispersion 0.58 Apixel -1 , R -1,800. Dichroic: 5,600 A, slit: 1".0. 



around Ha and the Nal doublet at 8184/8195 A. The 
Na I doublet clearly matches the velocity of the M-dwarf 
in the orbit, but the Ha emission seems to have a smaller 
velocity amplitude. A possible explanation for this is 
that the Ha emission comes from the M-dwarf surface 
that is closest to the WD, which may be heated by emis- 
sion from the white dwarf or the accretion region. 



Table 6 

Orbital Velocity Parameters of the M-dwarf. 



Parameter 


Value 


72 (kms -1 ) 


-48 ±5 


K 2 (kms -1 ) 


390 ±4 


to (HJD) 


54998.9375 ± 0.0003 


P (min) 


120.2090 ± 0.0013 



4.2. Photometric Variability 

The lightcurve of 1RXS J1730+03 shows clear period- 
icity (Figures [H - |3| , with two peaks per spectroscopic 
period. Most of the photometric data was acquired in the 
i' band, which contains contribution from both the M- 
dwarf and a cyclotron harmonic from the emission region 
near the WD. 

A Fourier transform of the data (Figure [T^|) shows a 
strong peak at sixty minutes. We interpret this as a 
harmonic of the orbital period. To determine the ex- 
act period, we analyze the data as follows. As the ob- 



ject has a short orbital period, we convert all times to 
Heliocentric Julian Date for analysis. We fit a sinusoid 
(mo + mASm([2n(t — t )]/P)) to each epoch, allowing mo 
and toa to vary independently for each epoch, but we use 
the same reference time to an d period P for the entire fit. 
The mean magnitude is correlated with the amplitude 
(Figure [5]). Note that the amplitude measured for the si- 
nusoidal approximation for each epoch is always less than 
the actual peak-to-peak variations of the source during 
that epoch, as expected. The source is in the active state 
in the first few epochs, and we exclude epochs 1-5 from 
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Figure 9. Blue component spectra of 1RXS J1730+03, showing the 12 late-time LRIS exposures (Section l2.2l Table IB") , The spectroscopic 
phase at mid-exposure is indicated for each spectrum. Dotted lines show spectra are obtained by subtracting a scaled spectrum of GL 694 
from spectra of this binary. The overlaid solid lines are smoothed versions of the same spectra. Vertical dashed cyclotron humps at 4440 A, 
5180 A, 6770 A & 8225 A. The cyclotron harmonic numbers (#3 — #6) are indicated in bold red. Two more cyclotron humps are seen in 
other spectra: a feature at 3540 A (Figure [6}, and a J- band feature (Figure 0. 



the fit, to avoid contamination from the accretion stream 
and/or the accretion shock. 

The best-fit solution is overplotted in blue in Figures 
[3] & [31 Since epochs 1-5 were not included in the fit, 
a sinusoid is overplotted in dashed green to indicate the 
expected phase of the variations. The best-fit period 
is 60.1059 ± 0.0005 minutes. This formally differs from 
the spectroscopically determined orbital period by 2.1a. 
However, this error estimate includes only statistical er- 
rors. There is some, difficult to determine, systematic 
error component in addition, so we do not claim any sig- 
nificant inconsistency. 

Periodic photometric variability for 1RXS J1730+03 
can be explained as a combination of two effects: cy- 
clotron emission from the accretion region and ellipsoidal 
modulation. The active state is characterized by a higher 
mass transfer rate from the donor to the WD, and results 
in higher cyclotron emission. This emission is beamed 
nearly perpendicular to the magnetic field lines, creating 
an emission fan beam. For high inclination systems, the 
observer crosses this fan beam twice, causing two high 
amplitude peaks per orbit (Figure [2]). In the active state 
when the emission is dominated by cyclotron radiation, 



the minimum leads the superior conjuntion by ~ 47°. 

In ellipsoidal modulation, and the photometric minima 
coincide with the superior conjunction. It is observed 
that in quiescence, the photometric minimum leads the 
superior conjunction of the white dwarf by ~ 14°. This 
suggests that the 0.29 ± 0.13 mag variation is caused by 
a combination of ellipsoidal modulation and cyclotron 
emission from the accretion region. 

4.3. Mass Ratio 

The semi- amplitude of the M-dwarf radial velocity, K2, 
gives a lower limit on the mass of the white dwarf. For 
a circular orbit, one can derive from Kepler's laws that: 



Mi 



PKl 
2ttG 



= 0.52 M, 



© 



(1) 



Tighter constraints can be placed on the individual 
component masse s Mi, Mg by co nsidering the geome- 
try of the system. lEggletonl (|1983[ ) expresses the volume 
radius i?2 of the secondary in terms of the mass ratio 
q = M 2 /M l : 



R2 

a 



0A9q 2 / 3 



0.6g 2 / 3 + ln(l + q 1 / 3 ) 



(2) 
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Figure 10. Velocity measurements 1RXS J1730+03. The best-fit 
solution gives P = 120.2090±0.0013 min, 72 = -48±5kms _1 and 
K 2 = 390 ± 4kms" 1 (Sectionl4~T) 



where the separation of the components is given by a = 
[P 2 G{M 1 +M 2 )/4tt 2 ] 1/3 • Thus for a fixed period P, R 2 
depends only on Mi. The radius of the white dwarf is 
much smaller than that of the M-dwarf. Hence, the M- 
dwarf will eclipse the accretion region if the inclination 
of the system is i < sin _1 (i?2/a). Figures El |4] show that 
we do not detect any eclipses in the system. 

Figuredlshows the allowed region for 1RXS J1730+03 
in a WD mass - M-dwarf mass phase space. The orange 
dotted region is excluded as the orbital velocity would be 
greater than the measured projected velocity. The red 
hatched region is excluded by non-detection of eclipses. 
The allowed mass of the primary ranges from the min- 
imum mass (Mi > 0.52 Mq) to the Chandrashekhar 
limit. The mass of the secondary is bounded above by 
the ZAMS mass for a M3 dwarf (M 2 = 0.38 M©). 

For a given Mi and a known orbital period, we can 
determine the radius of the secondary using Equation 
©, and can calculate the surface gravity (log g). For 
0.38 Mq > M 2 > 0.05 M Q , log g ranges from 5.1 to 4.8. 
This is consistent with log g — 5.0 for the best fit M- 
dwarf spectrum (Section l3.ip . 

The donor stars in CVs are expected to co-rotate. For 
1RXS J1730+03, the highest possible rotational velocity 
vsini is ~ 160 km s -1 for a 0.38 Mq, 0.27 R Q M-dwarf 
and a 1.1 Mq WD (Figure [13]). vs'mi will be lower if the 
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Figure 11. Velocity modulation of Ho and the 8184/8195 A Mai 
doublet. The spectra are offset by 2 X 10 -17 erg cm -2 s — 1 for clar- 
ity. The wavelengths are converted into velocities using the rest 
wavelengths 6563 A and 8195 A respectively. The vertical dotted 
line marks the radial velocity of the binary barycenter. The short 
dashed lines mark the radial velocities measured by fitting the com- 
plete spectrum. The Na I absorption follows the radial velocity of 
the M-dwarf, but Ha seems to have a smaller velocity amplitude. 
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Figure 12. Upper panel: Fourier transform of all i' data. Lower 
panel: expected Fourier transform for a pure sinusoidal variations 
with a period of 60.11 min, obtained by scaling and convolving the 
Fourier transform of the window function with a delta function 
corresponding to the best-fit period. 
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White Dwarf Mass (M ) 

Figure 13. The allowed range of masses (clear white region) for components of 1RXS J1730+03. The orange dotted region is excluded by 
the minimum inferred masses from radial velocity measurements. The red hashed region is excluded by non-detection of eclipses (Figures 
[3][4jl. The black dash-dot lines show contours of constant log g. The best fit spectra yeild log g = 5.0. The solid blue lines are calculated 
contours for the rotation velocity (kms -1 ) of the M-dwarf. A measurement of usini will help to constrain masses of the components. 



WD is heavier or if the M-dwarf is lighter. To measure 
rotational broadening in the spectra, we use a higher res- 
olution (R = 2 0, 000 ) template of the best-fit model from 
iZwitter et al.l (|2004l ). We take a template with zero ro- 
tation velocity and broaden it to differ e nt rot ational ve- 
locities using the prescription by Gra y (2005). Then we 
use our fitting procedure (Section [3.1|) to find the best-fit 
value for v sin i. For this measurement, we use only the 12 
relatively high resolution spectra from the second LRIS 
epoch (UT 2009 June 16). The weighted vsmi from 
the twelve spectra is 97 ± 22kms _1 , but the measure- 
ments show high scatter, with a standard deviation of 
54kms . We compared our broadened sp ectra with ro- 
tationally broadened spectra computed bv IZwitter et all 
(|2004D , and found that our methods systematically un- 
derestimate wsini by ~ 20kms~ 1 . We do not under- 
stand the reason for this discrepancy, hence do not feel 
confident enough to use this value in our analysis. A 
reliable measurement of v sin i will help better constrain 
the masses of the two components. 

4.4. Distance 

We estimate the distance to 1R XS J 1730+03 as fol- 
lows. Our fitting procedure (Section l3.ip corrects for ex- 
tinction and separates the WD and M-dwarf components 
of the spectra. We correct for varying sky conditions by 
using a reference star on the slit. The ratio of mea- 
sured flux to the flux of the best fit model atmosphere 



(T = 3500 K, log g = 5.0) is, 

/ 7 asurcd = ( = (6.1 ± 1.6) x 10- 23 (3) 

/model V " / 

where d is the distance to the source, and i?2 is given by 
Equation 

For 1RXS J1730+03, the maximum mass of the M- 
dwarf is 0.38 Mq, and the corresponding radius is i?2 = 
0.27 R®. Equation © then gives d = 800 ± 110 pc. This 
calculation assumes the largest possible M-dwarf radius, 
hence is an upper limit to distance. If the M-dwarf is 
lighter, say 0.1 Mq, we get R2 = 0.17 R®, yielding d — 
500±70pc. 

5. CONCLUSION 

1RXS J173006. 4+033813 is a polar cataclysmic vari- 
able, similar to known well-studied systems like BL Hyi, 
ST LMi and WW Hor in terms of the orbital period, mag- 
netic field and variability between active and quiescent 
states. This source is notable for the highly symmetric 
nature and high amplitude of the double-peaked varia- 
tion in the active state. This suggests a relatively high 
angle between the rotation and magnetic axes. Polari- 
metric observations of the source would help to better 
constrain the magnetic field geometry of the system. 

Most polars are discovered due to their highly variable 
X-ray flux. However, we mounted a followup campaign 
for 1RXS J1730+03 due to its unusual optical variability 
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properties. This suggests that current and future opti- 
cal synoptic surveys, such as PTfPI ()Law etall l2009l) 
and LSST can uncover a large sample of polars by cross- 
correlating opticaly variable objects with the ROSAT 
catalog. 
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